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ANALYTIC SHORT PERIOD LUNAR AND SOLAR 
PERTURBATIONS OF ARTIFICIAL SATELLITES 

David Fisher 

ABSTRACT 

The short period luni- solar theory of Kozai is gen- 
eralized for arbitrary obliquity of the ecliptic and 
inclination of the moon’s orbit to the ecliptic. Analy- 
tic first order lunar perturbations to the elements 
are derived. The theory is illustrated by an appli- 
cation to the communication satellite Intelsat 3F3. 
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ANALYTIC SHORT PERIOD LUNAR AND SOLAR 


PERTURBATIONS OF ARTIFICIAL SATELLITES 
1. INTRODUCTION 

In order to improve the accuracy of the computation of the orbits of satel- 
lites by analytic means, it is important to take into account lunar and solar 
perturbations . Thus the order of magnitude of the lunar perturbations of short 
period in the semi-major axis is 

S a/a = 1.1 x lCf 7 a 3 /a 3 (1) 

where a is the semi-major axis of the satellite and a e the radius of the Earth. For 
example, for a close Earth satellite with the semi-major axis equal to about 
7000 km, the order of magnitude of the perturbation will be about a meter, while 
for a communication satellite with semi-major axis of 42,000 kilometers, the 
order of magnitude of the lunar perturbations is about 1.4 kilometers. This is 
substantiated by the observations of an arc of Intelsat 3F3 data, Figure 1. 

In order to provide insight into our analysis, the luni-solar theories first 
developed by Kozai 1959, 1966 are briefly sketched. 

In 1959, Kozai showed that the luni-solar disturbing function R' is expres- 
sible in a form which we may write symbolically as 

R' = K(a, e, i) K' (a', e', i') cos (l, g, h, V, g' , h') (2) 

In Equation 2, the elements of the satellite are the unprimed quantities. Thus a, 
e, and i represent the semi-major axis, eccentricity, and inclination of the satel- 
lite orbit to the equator, while Z , g, and h refer to the mean anomaly, argument 



of perigee, and right ascension of the node of the satellite . The primed quanti- 
ties represent the corresponding elements of the sun or the moon. 

The elements i' and h*, when they refer to the inclination of the moon's 
orbit to the equator and the right ascension of the moon's node along the equator, 
show large non-linear variations with time. In 1966, Kozai completely elimin- 
ated i' in all terms and h' in most of the significant terms, by means of semi- 
analytic formulas of the type 

sin 2 i' cos 2h* = -156 + .065 cos N + .007 cos 2N (3) 

Here N is the longitude of the moon's node along the ecliptic. The quantity N is 
an angle which varies essentially linearly with time. Using formulas like those 
of Equation (3) the disturbing function is representable in the form 

R' = CK (a, e, i) cos (l, g, h, A', V , N) 

where C is a constant and A' is the mean longitude of the disturbing body. 

hi this report, the constants appearing in Equation 3 are replaced by more 
general formulas of the type 

s in 2 i' cos 2h’ = (e, J) + (e, J) cos N + (e, J) cos 2N (4) 

Here e equals the obliquity of the ecliptic and J the inclination of the moon's 
orbit to the ecliptic. Both e and J are very slowly varying quantities. 

Using formulas of this kind, the disturbing function of the moon is derived 
in a general form 
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R' = K (a, e, i) K' (a' , e' , e, J) cos (-£, g, h, V , N, V ) (5) 

The main advantages of using these transformations are that they can be 
applied to other three-body configurations than the Earth- satellite moon system, 
and that the accuracy of the transformation can be improved by using more pre- 
cise values of e and J. 

A program written in A PL language was prepared to evaluate the coefficient 
K’ of Eq. 5. The author used this language in a remote control unit for develop- 
ing the program, a listing of which is available at the Goddard Space Flight Center 
APL public library under the function name MNCLA. To my knowledge this is 
the first instance of a program in analytic celestial mechanics conducted ex- 
clusively from a remote control terminal. We learn to appreciate the remark- 
able capability of APL in handling vectors and matrices. 

Using this program for the coefficients K’, with the same initial conditions 
as Kozai it was found that Kozai’s development was correct with the exception 
of a single error in sign of a term of small magnitude. Thus the correct co- 
efficient of the term with argument fi- N - 2\' is +.003, rather than - .003 as 
given by Kozai. 

In this report, instead of applying the theory to close Earth satellites, we 
shall apply it to communication satellites, which are strongly influenced by the 
sun and the moon, to the extent indicated by Eq. 1. 

The results given in this report lead to the conclusion that the analytic 
formulas presented are accurate and convenient for obtaining the solar and lunar 
perturbations on a wide selection of satellites. 
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2. ANALYSIS 


The first step in the process of finding the analytic perturbations in the 
elements consists of finding and expanding the disturbing function. The next 
step is usually either to set up the Lagrangian equations of motion and solve 
these, or alternatively, we may define a determining function and use the method 
of canonical variables. In this report, the latter course is chosen. 

Following the procedure of Kozai, 1966, we first represent the dis- 
turbing function R® in the form 


/ J A vq A; q( c vq<q cos (q 0 + qV>* + ve ) (6) 

q.q' > v 

where q, q®, and v take the values given in Table 1, 

The nomenclature and symbols appearing in Eq, (6) are described in the 
section on notations, appendix 2. 

The functions A , A* . . and the constants c , are given in Table 1. The 
index q takes on the two values 0 and 2, the index q ! takes the value -2, 0, 2, 
while the index v varies over the range -2, -1, 0, 1, 2, as indicated in Table 1. 

The expansion of the disturbing function, Eq. (6) now begins. First, we give 
q® the values -2, 0, and 2 as indicated in Table 1. Then by means of Cayley’s 
tables, 1861, the disturbing function R' of Eq. (6) is expanded in terms of 
V , the mean anomaly of the disturbing body to obtain as in Fisher, 1971, 


R' = n' 2 m‘ a z ( - 
a 
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TABLE 1 


n 

■ 

H 


C Pq'q 

A '-V 

9A ^ 

B (cos i) 

0 

0 

0 

1 - 3/2 sin 2 i 

1/4 

1 - 3/2 sin 2 i’ 

3 cos i 

0 

2 

0 


3/8 

sin 2 i’ 


0 

0 


sin 2i 

3/16 

sin 2i* 

- 2 cos 2i/sin i 

0 

2 

1 , 


-3/8 

sin i 1 cos 2 i'/2 


0 

2 

I . 


3/8 

sin i' sin 2 i'/2 


0 

2 


sin 2 i 

3/8 

cos 4 i’/2 

- 2 cos i 

0 

0 

1 


3/16 

sin 2 i 1 


0 

2 

■i 


3/8 

sin 4 i'/2 


2 

-2 

0 

sin 2 i 

9/32 

sin 2 i' 

- 2 cos i 

2 

0 

0 


3/8 

1 - 3/2 sin 2 i' 


2 

2 

0 


9/32 

sin 2 i' 


2 

-2 


cos 4 i/2 

3/4 

cos 4 i'/2 

cos 2 i/2 

2 

0 



3/8 

sin 2 i' 


2 

2 



3/4 

sin 4 i'/2 


2 

2 


sin 4 i/2 

3/4 

cos 4 i'/2 

- sin 2 i/2 

2 

0 



3/8 

sin 2 i’ 


2 

-2 



3/4 

sin 4 i’/2 


2 

-2 

1 

sin i cos 2 i/2 

3/4 

sin i’ cos 2 i'/2 

1 - cos i - 2 cos 2 i 
2 sin i 

2 

0 

1 


-3/8 

sin 2i' 


2 

2 

1 


-3/4 

sin i' sin 2 i'/2 


2 

2 

-1 

sin i sin 2 i/2 

-3/4 

sin i' cos 2 i*/2 

(1 + cos i - 2 cos 2 i) 
2 sin i 

2 

0 

-X 


3/8 

sin 2i' 


2 

-2 

-1 


3/4 

sin i' sin 2 i'/2 



R* = n' 2 m* 



r 


< 


^O A v0 C vOO B O0 COS ^ 

+ A v0 A v0 c v00 B ['0 cos (i'^' +v9) 

+ A V0 A V2 C V20 B 1'2 COS + 2 g' + 

+ A v 2 A vO C v 02 B 00 COS ( 2 ^' + 

+ K 2 Ko C v02 B i'0 cos (i'^' + 2< £ + vB ) 

+ A v2 A v-2 C v- 22 B i'-2 cos (i'^' - 2 g' +24> + vd) 
+ K2K2 C v22 B i ' 2 COS C 1 *'*'* + 2 ®’ + 2t £ + 


(7) 


5 












The range of the indices i' and v are given in Tables 1 and 2. No confusion 
should arise between the use of i 1 as an index or i' as an inc lina tion. 


TABLE 2 


B \ » o 

D 

b\-_ 2 

D 

Ri 

^i'2 

D 

53/16 e’ 3 

-3 

845/48 e' 3 

-5 

1/48 e' 3 

B 

9/4 e* 2 

-2 

17/2 e* 2 

-4 

0 

0 

3/2 e* + 27/16 e’ 3 

-1 

7/2 e' - 123/16 e’ 3 

-3 

- 1/2 e ! + 1/16 e' 3 

1 

(1- e |2 r 3/2 

0 

1 - 5/2 e’ 2 

-2 

1 - 5/2 e’ 2 

2 

3/2 e s + 27/16 e’ 3 

1 

- 1/2 e' + 1/16 e* 3 

-1 

7/2 e' - 123/16 e’ 3 

3 

9/4 e’ 2 

2 

0 

0 

17/2 e’ 2 

4 

53/16 e’ 3 

3 

M 1 

> 

00 

1 

845/48 e’ 3 

5 


The coefficients B" lq> arise from the expansion by means of Cayley’s 
tables, Cayley, 1861, and are listed in Table 2 to order e' 3 . 

In order to simplify the notation we define 


Qo ^VO C v00 ®oo cos 


1*1 - c v02 B 00 cos (2^ + vd') 


+ A l 'o c ,oo B !o cos ( iT +^) +A' 0 c v02 B', 0 cos(iT + 20 + ^0) 

+ 2 g' + "*) + A ^,_2 c v _ 22 B'_ 2 cos - 2g’ + 20 + vQ') 

+ K 2 C v 22 B i2 COS + 2 S’ +20 + v6) 

Hence Eq» (7) is replaced by 

R ' = n' 2 m' a ! (l) 2 


( 8 ) 


( 9 ) 
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The indices and range of summation in this equation as well as those follow- 
ing are the same as in Eq. (7). 

By elementary trigonometry using the definition of 0, Eq. (9) can be ex- 
panded in the form 


R' 


n' 2 m* 



{ Ko Qo + A v 2 Qi cos 2 f - A„ a Qj sin 2 f} 


(10) 


Q 1 is obtained from P by replacing the variable <jfc by g, the argument of 
perigee. Q , is derived from Q x by replacing the cosine functions by sine 
functions, keeping the same arguments as in Q x> 

The quantities Qq , Q x , and Q x are no longer functions of the mean anomaly 
and are consequently slowly varying quantities which may be held constant when 
integrating with respect to the mean anomaly. 

To find the short period disturbing function, the mean value of R* found by 
averaging with respect to the mean anomaly is subtracted from R*. Consequently, 
the short period disturbing function R' is given’ by 



'2 


m 


V 

Z_, 




- A v2 Q, _ sin 2 f r , 




cos 2 f 


5 


e 


2 


( 11 ) 
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since the mean values of r 2 /a 2 , r 2 /a 2 cos 2f, and r 2 /a 2 sin 2f are 1 + 3/2 e 2 , 
5/2 e 2 , and 0 respectively. 

From the method of canonical variables applied to satellite theory,, 
Brouwer and Clemence, 1961, we define the determining function S by means of 
the partial differential equation 


1S_ 1 , 

n p 


(12) 


For arbitrary eccentricities, it is convenient to integrate Eq. (12) by transform- 
ing to the eccentric anomaly E. We use the relations 


r r 

— = 1 - e cos E, — cos f = cos E - e 
a a 


(13) 


dlrldE, 

a 


X* 


in f = (1 - e 2 ) 1/2 sin E 


Consequently, we find that 


with 


s = -'l-- J2 {A vo Q o T o ♦ o, T, + A v2 Q, T 2 } 


T = I - 2 e + - e 3 sin E + -e 2 sin 2 E sin 3 E 

0 V 4 / 4 12 


(14) 


T, = - - e 


sinE + i ^1 + ^-^sin 2 E - sin 3E 


(15) 


T 2 ~ 


(1 - e 2 ) 


2\l/2 


1 G 

cos E (1 + e 2 ) cos 2 E + — cos 3 E 

2 6 
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R ' of Eq. (11) can also be expressed in terms of eccentric anomaly since 

P 


r 2 o / 

cos 2 f =_e 2 - 2 e cos E + (1- — ]cos 2 E 

o2 2 l 2 


— sin 2 f = (1 - e 2 ) 1/2 (sin 2 E - 2 e) 

o2 


( 16 ) 


The perturbations in the Delaunay elements are then found by the formulas, 
Brouwer and Clemence, 1961, 


8 H = 


3 s 

8 l = 

3 S 

3^ 


3 L 

3 S 

K a - 

3 S 

3 g 

0 g - 

3 G 

3 S 

8 g = 

3 S 

3 h 


3 H 


(17) 


Substituting for S, using Eqs. (14) and (12), we have 


S L = - R ' 


2 n' 2 m' 


S G = L, {A - 2 Q 1 + A *2 T * } 


(18) 


H = CA -o Q ° T ° + A * 2 (Ql Tl + Ql T2>]} 


where Q 0 indicates the complement of Q 0 found by replacing the cosine functions 
in Q 0 by sine functions while preserving the arguments. 
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The perturbations in the angular elements are defined by 


8 ^ - 


8 g = - 


S h = - 


B S _ 

3 S 

B a 

B S B el 

B L 

[3 a 

B L 

Be B Lj 

B S _ 

/b s 

B e 

I- 

B S 

B G 

[3 e 

3 G 

B (cos i) 

B S _ 

B 

3 

B (co s i) 

B H 

B (cos i) 

B H 


3 (cos i) 

TcT 


(19) 


From the definition of S we have 


3 S _ 7 S 
B a 2 a 


(20) 


|S = £{a„ 0 Q 0 ♦ A„, Q, i h * A r2 Q, 


From the relations between the Delaunay and the Kepler elements we have 


B a/B L = 2/n a 

B e/B L = - [na 2 e(l - e 2 )] -1 

B e/B G = (n a 2 e) -1 (1 - e 2 ) -3/2 (21) 

B (cos i)/3 G = - cos i (n a 2 ) -1 (1 - e 2 ) _1/2 
B (cos i)/3 H = (n a 2 ) -1 (1 - e 2 ) -1/2 
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From the definitions of T , T , and T 2 , Eq. (15), we give the differential coef- 
ficients of these functions in the tabular form 


( 22 ) 


Using the functions in Table 1, and Eqs. (17)— (22) the perturbations in the 
elements may be calculated. It should be recalled that it is necessary to use 
mean elements of the satellite and the disturbing body as required by first order 
perturbation theory. For the case where either the sun or the moon is the dis- 
turbing body the appropriate elements are given in the Explanatory Supplement, 
1961. 

The element i’ for the sun is the obliquity of the ecliptic, a slowly varying 
quantity. The element h* for the sun is negligible. Consequently, the solar 
perturbations can be found by direct application of Eqs. (17) and (18). 
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However, the elements i' and h’ for the moon show large non-linear varia- 
tions. The analytic treatment for these elements is discussed in the next 
section. 

3. LUNAR PERTURBATIONS 

As shown by Kozai, 1966 the trigonometric relations useful in dealing with 
i ! and h' for the moon are given by 

cos i' = cos e cos J - sin e sin J sin N 

sin i' sin h' = sinj sin N (23) 

sin i' cos h' = sin e cos J + cos e sin J cos N 

Here i s and h’ are the inclination of the moon’s orbit to the equator and the 
longitude of the moon’s node along the equator respectively. The elements e , 

J and N are the obliquity of the ecliptic, the inclination of the moon’s orbit to the 
ecliptic, and the longitude of the moon’s node along the ecliptic. The elements e 
and J vary very slowly while N is essentially a linear function of time. 

As a consequence of Eq. (23), we find 

A; qf cos [fi - (q' + v) h']' 

(24) 

2 

= ( C vq'p + S ' S ',q'p) ( Sin V ) !P! COS (P + P N ) 

P=~2 
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where fi is an angle which does not contain h’. The integer s' is defined by 

s' = + 1 for (q' + v) p < 0 
s'=0 for (q'+r')p=0 
s ' = - 1 for (q' + v) p > 0 


For convenience we define the functions 


^vq'o ^vq'0’ 


S .q'o =S Wo =0 


for p / 0 C ' i = - C / 

^ vq p 2 v( l P 


S' - = 1 S , 

vq p 2 vc i P 


If we furthermore define the quantities d Q , d p and d 2 by the formulas 


1 

. 1 sin 2 e sin 2 T 

1 + 

2 (1 + cos e cos J) 

^ (l + cosecos J) 2 

sine 


2(1 + co se cos J) 2 


sin 2 -e 



4(1 + cos e cos J) 3 


we arrive at Table 3 



TABLE 3* 


C , and S , For Given Pairs v, q’ 

Vq'p Vq # p 7 ^ 


^vq*0 
^Vq' 1 
C «2 

(0,0) 

(1 - 3/2.sin 2 £ ) (1 - 3/2 sin 2 J) 
-3/2 sin 2 £ cos J 
3/4 sin 2 £ 

1 


*“Vq # 0 

<W1 

C|/q» 2 

(2,0), (0,-2), (0,2), (-2,0) 
sin 2 e (1 - 3/2 sin 2 J) 
sin 2 e cos J 
1/2 (1 + cos 2 £ ) 

«w. 

Svq'2 

2 sin 6 cos J 

COS € 

^Vq'O 

C , 

t'q t 
C *,' 2 

(1.0), (-1,0) 
sin 2 £ (1 - 3/2 sin 2 J) 
2 cos 2 e cos J 
- l/2 sin 2 € 

V. 

S Vq *2 

2 cos e cos J 
- sin e 

C ^q'0 

c ,, 

vq' 1 

(-1,2), (1,-2) 

1/2 sin e (cos J + cos e (1 - 3/2 sin 2 J)] * 
1/2 (cos 6 + cos 2 € cos J) 

- 1/8 sin 2 e 

Svq'l 

V. 

1/2 (1 + cos c cos J) 
- 1/4 sin e 

c , A 

vcf 0 

c ^> 

<W. 

c 

(1,2), (-1,-2) 

sin 3 £ cos J (1 - 5/2 sin 2 J) d 0 + 3/2 sin 2 £ cos £ (sin 5 J) d ; 
sin 3 £ cos J d ( + sin 2 £ cos £ ( 1+2 cos 2 £ ) d 0 
3/2 sin £ cos J (1 + cos 2 £ ) d 0 
+ sin 2 £ cos £ (1/2 + cos 2 J) d t + sin 3 £ cos J clj 

V. 

®vq' 2 

3 sin 3 £ d„ 

3/2 sin 2 t cos J d 0 + 3/2 sin 5 £ d. 

c 

vq' 0 

C r, 
Vq’ 1 

Cj/q'2 

(2,-2), (-2,2) 

1/4 (1 + cos e cos J) 2 + 1/8 sin 2 e sin 2 J 
- 1/2 sin e (1 + cos e cos J) 

1/8 sin 5 £ 



^Vq # 0 

C vq'l 

(2,2), (-2,-2) 

1/4 (1 - cos e cos J) 2 + 1/8 sin 2 € (1-32 cos 2 Jd 2 ) sin 2 J 
1/2 sin e (1 - cos € cos J) 

1/8 sin 2 € (1+32 cos 2 J d 2 ) 

Si/q'l 

^vq'2 

4 sin 3 € cos J dj 

4 sin 3 e cos J d 0 dj + l/2 sin 2 e cos € (1 + 2 cos 2 J) do 


*Th© subscripts vond q' aro chosen tp conform with Table 1. The numbers in parenthesis arc pairs of values of v and q . 


If we now apply Eqs„ (24) and (25) and use the values of C Vq(p and S Vq(p 
from Table 3 we may write R’ of Equation (9) in the form 

R* = n' 2 m' (l-j y ' A vq Q q » vi * p cos [q<£ + ^h + q'V + (i* - q') V + pN] (27) 
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where 


= c , B', , (C' , + s' S' , ) (sin J) 1 

Vq'q l'q' v Vq'p ^ yq ' p' K 


The coefficients Q q , vi , p » for the initial conditions e = 23!444 and J = 5°1454, 
were evaluated by the machine program discussed above, and are given in 
Appendix 1. 

The next step in obtaining the perturbations is to derive the determining 
function as outlined above. Corresponding to Eq. (14) we find for the moon 


s =. n,2 n m .. ' ..* 2 J2 {A w> Q q 'w q T o cos (^ h + + (i' - q‘) v + pN) 


[Tj cos (2oj + vh + q' V + (i' - q') V + pN) 


+ T 2 sin (2 co + vh + q'V + ( i ' - q') V + pN)]} 

T q , T , and T 2 and their differential coefficients are given by Eqs. (15) 
and (22) respectively. 

The lu nar perturbations of the Delaunay elements are then found by applying 
the canonical formulas given by Eqs. (18) and (19) to the lunar determining 
function S, given by Eq. (28). 

4. PERTURBATIONS IN THE SEMI- MAJOR AXIS OF INTELSAT 3F3 

A convenient way of illustrating the theory presented above is to con- 
sider the perturbations in the semi-major axis of Intelsat 3F3. Osculating 
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elements of this satellite were available over a 21-day arc, nearly all of the 
observations being taken at 2-day intervals. 

The definition of the Delaunay element L is given by 

L = (jjl a) 1/2 = na 2 (29) 


From Eq„ (29) we find 


<- 2na 2 „ T 

Sa = _SL (30) 

Using Eq. (18), the Vanguard system of 'units where fj. = 1, and the Earth's 
equatorial radius a e as a normalizing factor, we have 

Sa/a e = 2 R; (31) 

The eccentricity and inclination of Intelsat 3F3 are about .005 and 0°34 s 
respectively. This permits us to derive a very abbreviated expression for the 
short period disturbing function which is adquate for the semi-major axis. Thus 
we find using Eq. (27) and the results of appendix 1 that 


Sa/a = 1.102 x 10" 7 J cos 4 i (.682 cos {2l + 2g + 2h - 2\') 
+ .131 cos {2l + 2g _ 2\’ - V)} 


(32) 
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The factor 1.102 x 10~ 7 arises from 2n ,2 m' for the moon in Vanguard 
units. 

The perturbations obtained by applying Eq. (32) to the Intelsat 3F3 data are 
plotted in Figure 1. Agreement between the observed and calculated values of 
the semi-major axis is adequate, since only the major perturbations are 
considered. 

It should be recalled that Intelsat 3F3 is a synchronous satellite with the 
period of a day. Since the satellite was observed at equal multiples of a day, 
the mean anomaly i increased by multiples of 277 for each observation, thus 
appearing as a constant angle in Eq. (32). Consequently as shown by Eq. (32), 
the perturbations in the semi-major axis should show a dominant semi-monthly 
period. This is borne out by the observations, (see Fig. 1). 

i 

5. RESULTS AND CONCLUSIONS 

A general first order short period theory is given for obtaining perturba- 
tions analytically either for the earth, satellite, moon, sun system or a configura- 
tion of a similar type. 

A computer program was written in APL to apply this theory. The program 
was used to compute coefficients needed for the disturbing function for the spe- 
cific case of e = 232444 and J = 521454, and the results shown to be in general 
agreement with the work of Kozai, 1966. 

The theory is illustrated by computing the perturbations in the semi-major 
axis of Intelsat 3F3 with good agreement obtained with observations. 
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The analysis shows that short period lunar perturbations vary from about 
one meter for close earth satellites to about 1.4 kilometers for synchronous 
type satellites. The lunar perturbations as well as those due to the sun can be 
computed analytically by the formulas of this report. 
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SEfVil-MAJOR AXIS (KILOMETERS) 



Figure 1. Semi-major axis vs time for satellite intelsat 3F3- 
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APPENDIX 2 


NOTATIONS 

Ay function of i 

A^, q , function of i’ 

A' , function of e and J 

vq p 

B. function of e 

B' , , function of e’ 

1 q 

C Vq , function of e, J 

G Delaunay variable = (^ a (1 - e 2 )) 1/2 
H Delaunay variable = (/xa (1 - e 2 )) 1/2 cos i 
j Inclination of moon’s orbit to ecliptic 
L Delaunay variable = (yua) 1/2 
N Longitude of lunar node along ecliptic 
R' disturbing function 

S < function of e and J 

a semi-major axis of satellite 
a e mean equatorial radius of earth 
c y q » q constants of the disturbing function 
e eccentricity of the satellite 
e* eccentricity of the disturbing body 
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f true anomaly of satellite 
g argument of perigee of satellite 
g' argument of perigee of disturbing body 
h right ascension of node of the satellite 
h' right ascension of the disturbing body 
i inclination of satellite orbit 
i' inclination of the disturbing body 
l mean anomaly of satellite 
V mean anomaly of the disturbing body 

m ' ratio of mass of the disturbing body to mass of earth plus mass 
of the disturbing body 
n mean motion of satellite 
n' mean motion of disturbing body 
p summation index 
q summation index 
q' summation index 
r geocentric distance of satellite 
r' geocentric distance of the disturbing body 
e obliquity of ecliptic 
4) argument of latitude of satellite = f + g 
<56' argument of latitude of moon = f’ + g’ 

K' mean longitude of the disturbing body = -V + g* + h' 
fj. gravitational constant 
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